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Abstract 

A systematic gauge-invariant method, which starts directly from QCD, 
is used to calculate the rate for an upsilon meson to decay inclusively into 
a prompt photon. An expansion is made in the quark relative velocity f , 
which is a small natural parameter for heavy quark systems. Inclusion of 
these O(w^) corrections tends to increase the photon rate in the middle z 
range and to lower it for larger z, a feature supported by the data. 



INTRODUCTION 



The hadronic decays of the T family of bb mesons proceeds mainly via an inter- 
mediate state consisting of three gluons. By replacing one of the outgoing gluons 
with a photon one obtains the leading order contribution to the production of 
direct photons, i.e. those photons which do not result from 7r° decay, etc. The 
spectrum of such photons provides, in principle, an excellent test of perturbative 
quantum chromodynamics (QCD) because in this case one has a large number 
of data points against which theoretical predictions can be compared. This is 
in contrast to the prediction of a decay rate, which is a single number. How- 
ever it is well knownlH that the photon spectrum ^Si — > 7 + X, calculated at 
leading order |0, is too hard when compared against experiment, both in J/\E' 
and T decays. Such calculations yield an almost linearly rising spectrum in 
z = 2E^/M with a sudden decrease at z = 1. A next-to- leading order calculation 
by Photiadis[3 sums leading logs of the type ln(l — z) and yields some softening. 
However, the peak is still too sharp and close to z = 1. An earlier calculation by 
Field 1^ predicts a much softer spectrum which fits the relatively recent data^] 
quite well. This uses a parton-shower Monte Carlo approximation wherein the 
two gluons recoiling against the direct photon acquire a non-zero invariant mass 
by radiating further bremsstrahlung gluons. This does not, therefore, qualify it 
as an ab-initio perturbative QCD calculation. We note that in refs[2-4] the non- 
perturbative dynamics of the decaying hadron is described by a single parameter 
0(0), the quark wavefunction at the origin. This leads to the assertion that the 
ratio of widths for decay into prompt photons and l^l^ pairs is independent of 
quark dynamics. 

In this paper we compute the rate for ^Si ^ 'j + X taking into account 
the bound state structure of the decaying quarkonium state. We note that the 
description of hadron dynamics in this decay process by just 0(0) is correct only if 
one assumes that Q and Q are exactly on-shell and at rest relative to each other. 
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This assumption is only approximately true - heavy quarkonia are weakly bound 
QQ composites and f^/c^ is a small parameter. Improvement requires introduc- 
tion of additional hadronic quantities, which we identify within the context of a 
systematically improvable gauge-invariant theory for quarkonium decays. This 
formalism has been recently applied to one and two particle decays^, 0. Here 
we apply the method of ref0 to the more complicated three particle case and 
obtain the photon spectrum for the process T — 7 + 2g. We find that inclusion 
of binding and relativistic effects via the two additional parameters, e^/M and 
V^0(O)/M^0(O), makes the computed spectrum softer for large z, {z < 0.9). For 
still larger z, 0.9 < z < 1, there are non-perturbative effects due to final-state 
gluon interactions which cannot be reliably computed and which, therefore, we 
shall not address. 



FORMALISM 

The starting point of our approach is that the decay amplitude for ^5*1 7 + X is 
given by the sum of all distinct Feynman diagrams leading from the initial to the 
final state. The first step is to write a given diagram in the form of a (multiple) 
loop integral. Consider, for example, one of the six leading order diagrams (Fig. 
la). Omitting colour matrices and coupling constants for brevity, its contribution 
can be expressed as 

M[k) is the usual, but obviously non-gauge invariant zero-gluon, Bethe-Salpeter 
amplitude, 

M{k) = I d^xe"'-''{0\T[^P{-x/2)ij{x/2)]\P,e). (2) 

In equations 1-2, is the relative distance between quarks, fc^ is the relative 
momentum, = M^, and Si = qt — \P. We define the binding energy as 
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eB = 2m — M. Provided all propagators are far off-shell, they may be expanded 
in the two small quantities es/M and k/M. This yields the expression, 

+ (OlV'^Ha idp ^\P, e)^9"9'3/i^i^2^3 + ...]. (3) 

We have defined d = — d ), and /i^i/^^Ms jg "hard part" which combines 
terms from all six leading diagrams^. One can readily see that it is the sum of 
terms of the type in eq.l corresponding to different permutations of indices and 
momenta. There are 12 one-gluon diagrams one of which is illustrated in fig lb, 
which must be added as corrections to the no-gluon amplitude. These all have 
the general form 

Tr^'' = J (04^^^ ^''(^' k')H^^^^^^{k, k'), (4) 

where Mp{k, k') is a generalized B-S amplitude, 

Mf{k,k') = J d1rd'ze*'=-V'='-"(0|T[V'(-x/2)A^(z)^(x/2)]|P,e). (5) 

The gluon which originates from the blob is part of the QQg Fock-state compo- 
nent of the meson and has its momentum k' bounded by R^^ ^ k' <^ M, where R 
is the meson's spatial size. Hence it is to be considered soft on the scale of quark 
mass. Again, one may expand the propagators in Hj^^'^^^^{k, k') about k = k' = 
to get 

j,mM2/.3 ^ TrlMf'H^'''^^''' + MP^'^d^H^'^'^^''-' + M'^^'^d'^H^'^'^'''' + ...], (6) 

where, 

Mf = {OltpipAPlP^e) 

M'^'° = {0\^^idAP\P,e), (7) 
^ In actual fact, only three of the diagrams need to be evaluated because of time-reversal 
symmetry. This simplification halves the number of diagrams in the one-gluon and two-gluon 
cases as well. 
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The derivative id acts only upon the quark operators. 

The two soft-gluon contribution (figs.lc-ld) to the amphtude is handled sim- 
ilarly but is more complicated. The matrix elements which enter into T2^^'^^''^{k) 
are of the type < 0\iljiljAPA^\P,e > and so on; we shall not list these explicitly 
here since the principle is rather clear. It is straightforward to show that in the 
sum Tq + Ti + T2, all d's combine with A's to yield covariant derivatives and/or 
field strength tensors, 



(To + Ti + Ta)^^'^^^^ = Tr[(0|V^|P, e)/^^^'^^^^ + {0\ipiDa'^\P, e)d''h''' 
+ {0\4!iDa^Dp^p\P,e)^^"^'^h'''^'^^'^ + {0\^F"^^\P,e)^d'^H^'''"'' + ...]. (8) 

The above is a sum of terms, each of which is the product of a soft hadronic 
matrix element and a hard perturbative part. 

To proceed, one can perform a Lorentz and CPT invariant decomposition of 
each of the hadronic matrix elements in Eq 10. This is somewhat complicated |^ 
and involves a large number of constants which characterize the hadron. Consid- 
erable simplification results from choosing the Coulomb gauge, together with the 
counting rules of Lepage et. al [|]. The upshot of using this analysis is that, in 
this particular gauge, the gluons contribute at 0{v^) to the reaction '-f + X, 

and hence can be ignored. Even this leaves us with too many parameters, and 
forces us to search for a dynamical theory for the 1 quarkonium state. We shall 
assume, in common with many other authors, that the Bethe-Salpeter equation 
with an instantaneous kernel does provide an adequate description. This has been 
conveniently reviewed by Keung and Muzinich|]TO[ and we adopt their notation^. 
The momentum space B-S amplitude x{p) satisfies the homogeneous equation, 

Xip) = tGoiP,p) J ^ K{P,p,p') x{p), (9) 



^We find the analysis of rcf to be wanting because it does not properly deal with the issue 
of gauge-invariance of the meson state. Further, while the binding energy is taken into account, 
the wavefunction corrections - which are essentially short-distance or relativistic effects - are 
not. 
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which, after making the instantaneous approximation K{P,p,p') = V{p,p') and 
reduction to the non-relativistic hmit yields, 

Mi/2(M -2E){E + m- p. 7) /(I - 7o)(^ + m - p.^)^{\p\) 



xip) 



AE{E + m){p^ + f - E + ie){p° - f + E - le) 



(10) 



The scalar wavefunction 0(|p|) is normalized to unity, 

d^p 



and. 



E 



fp + m?. 



:iii 



(12) 



Fourier transforming x{p) to position space yields (0|'(/'(— a;/2)^/'(a;/2)|P) from 
which, by tracing with appropriate gamma matrices, the coefficients below can 
be extracted. So finally, to O(f^), one has a rather simple result. 



(0|^^|P,e) 
(0|V'ia„Vl^,e) 



^mV^ f 1+4 1+^1 /-V/^^ 

2 \ M^I^\ mI^ 2 3M2 



3 M2 
6 M2 



pi/ 



M2 



(13) 



DECAY RATE 



All the ingredients are now in place to calculate the decay T ^ 7 + 2g. In 
squaring the amplitude obtained by substituting Eqs.l3 into Eq.8, terms involving 
the product of and may be neglected. We assume the emitted gluons to 
be massless and transverse, and to decay with unit probability into hadrons. 
Polarizations of the final-state particles are summed over, and the spin states of 
the initial meson are averaged over. Summing over final-state colors yields 2/3, 
and one must include a factor of 1/2 for identical gluons. The Lorentz invariant 



phase-space factor for 3 massless particles has a standard expression] 13] which is 
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best expressed in terms of the dimensionless energy fractions Xi = 2Ei/M which 
satisfy xi + X2 + = 2. The variables s, t, u are symmetric functions of Xj, 



M2(1-xi), 



(14) 



Ignoring radiative radiative corrections for the moment, a tedious calculation^ 
yields, 



256 



dxi dx2 9 



2 2 I 



[^o/o(s, u) + iiEfsis, t, u) + riwfw{s, t, u)] . (15) 



e„ is the quark charge and. 



The function /q provides the standard, leading order result: 

(S^ ^2 + ^2 ^2 _^ ^2 g2 ^ ^2 _g ^ ^\ 



(16) 



(17) 



(S-M2)2(t-M2)2(u-M2)2 ' 

The binding energy and wavefunction corrections, fs and fw respectively, are 
more complicated: 



+ (s^t^ + t^u^ + u^s^) (s^ + + + 15 s t 
+ stu[s^ + + u^) +29 M^sH\'^ 



u 



fw{s,t,u) 



Ulstu + t^ + M^) - 85 (^sh^ + t^u^ + u^s^ 



27 ( sH'^ + t^u^ + M^s^ ][s'' + t'' + u 



205 

2 „2j.2„,2 



stu 



139 s t M (s^ + + M^) - 463 52^2^ 



:i8) 



The denominator D is, 



= (s - M2)^(t - M''f{u - M''f. 



2\3/ 



r2\3 



(19) 

We used Mathematica||T^ , supplemented by the HIP package [p^, for computation of traces 
and simplification of algebra 
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Integrating over the energies of the outgoing gluons for a fixed photon energy 
yields 

dV 256 , o|0(O)P , ^ , , ^ , , ^ . M 

— = —e^a.a^^^ ^^^p^^^^ ^ r^BFB{z) + r^wFw{z)] , (20) 

where z = 2E^/M and, 

Fo = [l+4^-2e3-e^-2e5 + 2e(l + 2e + 5aiog^]/(l-0'(l + 0', 

= [2 - 16^ + 10^2 - 48^=^ - 10^^ + 64^^ - 2^^ 

+ (1 - 3^ + 14^' - 106^^ + 17^^ - 51^^) log^]/2 (1 - + 0', 

Fw = [-26 + 14^ - 210^2 ^ 134^3 ^ 274^^ - 150^^ - 38^^ + 2 xi'^ 
- (27 + 50^ + 257^2 - 292^3 ^ 205^^ - 78^^ - 41^^) log^]/ 

3(l-0'(l + 0'- (21) 

In the above, ^ = 1 — z. The integrated decay width isQ, 

128, n , o o|</)(0)P f as \ . . 

= — (vr' - 9)ejaea^^^^ [l + a-^ - l.OSriB + l9Mr]wj • (22) 

Where we have included the radiative corrections of 0{as) which are of the same 
order in f^/c^ as the other corrections, and were calculatedO] many years ago, 



a = /3oln(/i/mQ) -4.37- 0.77n/, (23) 

where jSo = 11 — 2r2j/3. The parameters rjw and tjb are independent of each 
other in the present treatment. We note, however, that if we impose the condi- 
tion r]w = ^TjB then the result Eq. 3.5 of Keung and Muzinich[10] is precisely 
recovered. This latter condition is equivalent to -gV^0(O) = |eB0(O), which is 
the Schrodinger equation for quark relative motion in a potential which vanishes 
^Note that Eq.22 does not take into account non-pcrturbative effects citeVoloshin which are 
significant in the part of the phase space where one of the quark propagators become soft, and 
where the gluon vacuum condensate plays a role. 
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at zero relative separation. It is also worthy of note that the same condition 
emerges as a renormalization condition in the treatment of positronium by La- 
belle et al|T5[ (see their equations 11 and 12). However in our treatment there 
is no principle which apriori constrains r]B to bear a fixed relation to rjw and 
therefore both will be considered adjustable parameters. 

The application of Eq.22 must be done with caution because extraction 
of the direct photon decay rate from the data requires an extrapolation down 
to small photon energies. But in this energy range the prompt photons are 
heavily contaminated by photons from 7r° decays. A numerical estimate of the 
correction factors requires the value of tib and ri\Y. We have chosen rrib = 4.5 
which gives t]b = —0.048. If we take as = 0.20 then rjw can be fixed by using 
the experimentally known numbers|17[j, 



r(T^2^ + 7) = 1.28 ±0.10 if eV. 

r(T^/r) = 1.34 ± 0.04 ifeV. (24) 

This gives a range of values for rjw- We have plotted the graphs in fig. 2 at rjw = 
—0.0059. The binding, Fb{z), and wave-function, Fw^z), correction terms tend 
to cancel each other over part of the z region. The effect of final-state interaction 
corrections can be reasonably well estimated [3] provided one stays away from 
the end-point z = 1. In fig. 2 we compare the data, taken from Ref[5], with the 
prediction of our model appropriately folded with the experimental photon energy 
resolution (assumed to be Gaussian). The effect of the binding and wavefunction 
corrections calculated in this work is sizeable, and tends to to increase the photon 
rate in the middle z range and to lower it for larger z. While this appears to be 
in the right direction, it would be highly desirable to have more precise data. 
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SUMMARY 



The approach taken in this paper for calculating the amplitude for T — > 7 + X is 
to take the sum of all distinct Feynman diagrams leading from the initial quarko- 
nium state to the final state. Each diagram is put into the form of a (multiple) 
loop integral with a kernel which is a product of a hard part and a soft part. The 
hard part is treated with perturbative QCD, and the soft part is analyzed into 
its different components with the use of Lorentz, C, and V symmetries. Use of 
the QCD equations of motion enables separation of these components according 
to their importance in powers of v. At the last step, a specific commitment to 
dynamics is made and the B-S equation is used to express the components in the 
form of wavefunctions. However, the un-regularized value of V^0(O) is singular 
at the origin V^0(O) ~ M(j){0)/r. As is clear from the uncertainty principle, 
the local kinetic energy becomes very large at short distances and the expan- 
sion in powers of v breaks down. This difficulty was circumvented by imagining 
that annihilation takes place in a diffused region of size 0(l/m), i.e., that 0(0) 
and V^0(O) are quantities renormalized at this scale, and to be considered as 
adjustable parameters. The numerical investigation we undertook showed that 
varying these within reasonable limits led to substantial improvement in the in- 
termediate z region but was insufficient to reproduce the data near z — 1, once 
again underscoring the importance of final-state interactions between coUinear 
gluons. 
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Figure Captions 



1. a) One of the six leading order diagrams. 

b) One of the 12 one-gluon diagrams. 

c) One of the 24 two-gluon diagrams . 

d) One of the 12 gluon self-couphng diagrams. 

2. The photon spectrum as a function of folded with the experimental photon 
energy resolution. The dotted line is the zeroth order QCD result, the dashed 
line incorporates the binding and wavefunction corrections, with r]B — —.048 and 
Tjw = +0.0059 . The solid line is the final result including final-state interaction 
of Ref[3]. 
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